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ABSTRACT 

We have investigated the crustal properties of neutron stars without fallback accre- 
tion. We have calculated the chemical evolution of the neutron star crust in three 
different cases (a modified Urea process without the thermal influence of a crust, a 
thick crust, and a direct Urea process with a thin crust) in order to determine the 
detailed composition of the envelope and atmosphere as the nuclear reactions freeze 
out. Using a nuclear reaction network up to technetium, we calculate the distribution 
of nuclei at various depths of the neutron star. The nuclear reactions quench when the 
cooling timescale is shorter than the inverse of the reaction rate. Trace light elements 
among the calculated isotopes may have enough time to float to the surface before 
the layer crystallizes and form the atmosphere or envelope of the neutron star. The 
composition of the neutron-star envelope determines the total photon flux from the 
surface, and the composition of the atmosphere determines the emergent spectrum. 
Our calculations using each of the three cooling models indicate that without accre- 
tion of fallback the neutron star atmospheres are dependent on the assumed cooling 
process of the neutron star. Each of the cooling methods have different elements com- 
posing the atmosphere: for the modified Urea process the atmosphere is 28 Si, the thick 
crust has an atmosphere of 50 Cr, and the thin crust has an atmosphere of 40 Ca. In all 
three cases the atmospheres are composed of elements which are lighter then iron. 
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Neutron stars are the end product o f the core collapse of 
a sta r with a mass greater than 8Mia IjLattimer fc Prakashl . 
|2004 ). A canonical neutron star has a radius of 10 km and a 
mass of 1.4Mq, resulting in an average density greater than 
that found in an atomic nucleus. A neutron star is com- 
prised of five major region s: core, inner crust, outer crust, 
envelope, and atmosphere l|Lattimer fc Prakashl . |2003 ). The 
crustal layer extends 1-2 km below the surface inward to a 
density of around 10 14 g/cm 3 . Here we define the neutron- 
star envelope to be the upper layer of the crust that throt- 
tles the heat flow running from a density of 10 7 g/cm 3 
outw ard (|Hernquist fc Applegatel . 1 1984 iHevl fc Hernquistl 
l200ll ). The atmosphere lies at relatively low density and 
comprises a column density of about 1 g/cm 3 . 

The envelope and atmosphere contain a negligible 
amount of the neutron-star mass (or the mass of the crust 
for that matter) but play a crucial role in shaping obser- 
vations of neutron stars. In order to fully interpret the ob- 
served emission, we need to understand the crustal compo- 
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sition especially the lightest trace elements that can float 
to the surface to form the envelope and atmosphere. The 
atmosphere shapes the emergent spectrum, understanding 
the composition could yield predictions of spectral features 
in the neutron star thermal emission. The envelope, in turn, 
influences the transport and release of thermal energy, a 
change in its composition changes the thermal conductiv- 
ity and inferred surface te mperature of the neutron star 
IlLattimer fc Prakashl [2001 ). 

Without fallback from the supernova or other accre- 
tion, a neutron star is expected to have a n atmosphere of 
iron- group elements l|Chiu fc Salpeterl , Il96 4) . An example of 
a family of neutron stars which are closest to this ideal- 
ized situation are the isolated neutron stars. The isolated 
neutron star RX J185635-3754 has been studied extensively. 
The surface composition of RX J 185635-3754 has been ex- 
amined by fitting the x-ray spectra to various model atmo- 
spheres, including models attributed to fallback accretion. 
These models have included: black body, hydrog en, helium, 
iron, silicon-ash atmospheres l|Pons et all I2002T ). and later 
extended to two black bodies, pure silicon, low iron sili- 
con a sh and magnetic hydrogen atmospheres l|Walter et all 
2004). The model atmospheres that fit the spectrum 
of RX J 185635-3754 the best are those with heavy el- 
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eme nts, though the predicted absorptio n lines are not 
seen (|Pons et all |2002| ; Iwalter et all . |2004 ) . 

Here we calculate the mass fractions that exist 
in the crust after a neutron star cools long enough 
for the nuclear reactions to be quenched. These 
mass fractions are calculated using a 489 isotope 
reaction network, torch (the code is available at: 
http : / / cococubed.asu.edu/ code_pages/net_torch. shtml 
iTimmesl . |l_999j)- The results from the modified Urea case 
are compared to analytic calculations of the freeze-out 
of the neutron star crust using the nuclear statisti- 
cal equilibrium software nse (the code is available at: 
http : / / cococubed . asu . edu/ code_pages/nse . shtml ) and 
the reaction rates in torch. 



2 NEUTRON STAR COOLING 

In order to calculate the crustal abundances three different 
cooling curves were used: a modified Urea process, a thick 
crust, and a thin crust. Each of the cooling curves start near 
a temperature of 10 10 K and cools until the nuclear reactions 
are quenched. 

In the case of the modified Urea process the nuclear 
reactions are quenched before a year has past. The core 
temperature and the c ooling time are determined by the 
modified Urea equationfSha piro fe Teukolskvl Il986l ): 



At = lyrT 9 - 6 (/) { 1 



21.(0 



(1) 



where Tg (/) is the temperature of the outer core (T c ) in units 
of 10 9 K. For simplicity we have assumed that the tempera- 



g/cm 



ture of the crust is isotherma l at densities above 10 7 
(e.g. iHevl fc Hernquistl 1200 ll ). 

As the crust of the neutron star cools via the diffusion 
of heat the to interior, the surface of the star is expected to 
remain ho t until the hea t reser voir of the crust has been 
exhausted l|Lattimer et all I1994T ). The thick and the thin 
crust cooling models take into account this thermal discon- 
nect between the crust and core, whereas the modified Urea 
model does not. 

The cooling curve f o r a t hick crust is based on the 
model in lLattimer et alj (| 19941 ). In this model the core is 
cooling rapidly and the cooling wave takes 15 years to dif- 
fuse through the crust. This model is appropriate for a nor- 
mal neutron star. The nuclear reactions are quenched be- 
fore the the cooling wave passes through the crust, thus the 
core cooling does not affect the reactions in the crust. The 
cooling curve used in this work is taken from Figure 3 in 
lLattimer et al.l (|l994l ). where the hottest parts of each curve 
set the age and temperature of the neutron star. 

The thin crust cooling curve is appropriate for a 
strange star. In this case the core is cooling via the 
direct Urea pr ocess, where the tem perature relation is: 
t ~ 20 r c T 4 s (|Lattimer et all . Il99ll ). The crust of the 
neutron star is cooled by crustal bremsstrahlung. The 
cooling of the crust is determined from the crustal lu- 
minosity and the total thermal energy of the crust. 
The luminosity of crustal bremsstrahlung is: Li rems ~ 
(5 x 10 39 erg/s)(Af cr /M Q ), where M cr is the mass of the 
neutron star crust and the total thermal energy of the 



crust is: U cr = \kT^- (|Shapiro fc Teukolskvl Il986l ). 
where A is the average atomic mass and m u is the mass 
unit. The crust of a strange star reaches a density of a 
few times 10 10 g/cm 3 , ins tead of neutron dr i p in t he case 
of a normal neutron star l|Steiner fc Madsenl . l2005h . For a 
crustal densit y of 4.8 x 10 10 g/cm 3 the equilibrium nu- 
cl eus is 80 Zn(fshapiro fc Teukolskvl 1 19861 ). From Figure 1 
of lBrown et alT ( 19981 ) the time scale for the cooling wave to 
pass through the crust is on the order of a month. As a re- 
sult, for the thin crust case the crustal cooling is dominated 
by the crustal bremsstrahlung for the first month. After a 
month has past the direct Urea process dominates. We take 
the temperature of the crust to be 



T — 2brems (e 



-(t/month)' 1 



T DU (1 



-ft/month)"* 



(2) 



where Ti, rems is the temperature of crustal bremsstrahlung 
and Tdu is the temperature at a specific time step us- 
ing the equation for direct Urea cooling. The T(, rems starts 
at a temperature of 10 10 K and decreases by 0.1% each 
step. The time, or age, of the star is given by the crustal 
bremsstrahlung equation for a thin crust: 



t = J^L 2 x 10- 31 ^T 9 - 5 [1 - 2f l 
1 ' ' erg 



10 Amn 



(3) 



for a neutron starstarting at a temperature of 10 10 K, where 
T10 is the temperature in units of 10 10 K. 

The cooling curves for the 3 cases are shown in Fig- 
ure [1] We have assumed for densities above 10 7 g/cm 3 the 
crust is isothermal and Figure [1] reflects the cooling curves 
appropriate for densities at 10 7 g/cm 3 and greater. For den- 
sities below 10 7 g/cm 3 the input temperature is interpolated 
between the core (T c ) and the surface (T a ): 



log(T) 



log(Te) - log(T a ) 



7 



Iog(p) + log(2V 



(4) 



where the su rface temperature i s giv en as: T s — 
(10 K 1/2 T C ) 2/3 (|Shapiro fc Teukolskvl mm . 



3 CRUSTAL MASS FRACTIONS 
3.1 Running torch 

The code torch is a reaction network written by 
Frank X. Timmes and is available on his website 



(http://cococubed.asu.edu). The software follows the 
abundances of 489 isotopes to technetium. T he construction 
of the reaction network is covered in detail in lTimmes et al.l 
(2000) and summarized below. The reaction network starts 
by determining the mass fraction of an isotope, i,: Xi — pi/p 
and the corresponding molar abundances of the isotope: 
Yi = Xi/Ai. A set of partial differential equations are con- 
structed from the continuity equation of the isotope: 



HY 

^ + V ■ (YM) = Rr 



(•>) 



where R\ is the total reaction rate and Vi is the mass dif- 
fusion velocity, which is set to zero. This results in a set 
of ordinary differential equations that comprise the reaction 
network: 



dY _ £ 
~dt ~ 



(6) 
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Figure 1. The three different cooling curves used for the abun- 
dance calculations. The methods for cooling include the modi- 
fied Urea, thick crust (labelled Lattimer 1994) and a thin crust 
(labelled Strange Star). These curves are appropriate for densi- 
ties at 10 7 g/cm 3 and above. For densities below 10 7 g/cm 3 , the 
curves would be interpolated between the surface temperature 
and the temperature at 10 7 g/cm 3 . The thick crust is taken from 
lLattimer et al.l £1994), and the thin crust is appropriate for a 
strange star. 



The system of ordinary differential equations is integrated 
using a Bader-Deuflhard method coupled with MA28 sparse 
matrices. Three methods for inte grating and eig ht different 
matrix packages are compared in iTimmea (|l99S ). 



3.2 Mass Fraction Calculations 

For each of the densities investigated the system is started 
in nuclear statistical equilibrium. The temperature and the 
length of time at which the system burns is determined by 
the cooling curve. After the reactions have quenched the 
abundances at a specific density are calculated and the iso- 
topes which are sufficiently abundant to compose the atmo- 
sphere or envelope are determined. 

In order for the isotope to form the atmosphere it needs 
to have a surface density greater than lg/cm 2 ~ ar/nip, 
the ratio of the Thompson cross-section to the mass of the 
proton. The first step in determining which isotopes rise to 
the surface is to calculate the total pressure at the density 
of the calculation. For densities greater then 10 6 g/cm 3 the 
pressure is given by: 



P = 



1.2435 x TO 15 



4/3 
Me 



lg/cm" 



4/3 



dyne/cm 



(7) 



where we have assumed that relativistic electrons dominate 
the pressure and that fi e = 2. For the case where we look at 
a density of 10 6 g/cm 3 the pressure is given by: 

2 



4>(x) = [x^l + x 2 {2x 2 /3 - 1) + ln(a; + y/l + x 2 )] 



(9) 



and x = pf/m e c, the ratio of Fermi momentum (pf) to 
the product of the mass of the electron and the speed of 
light (|Shapiro fc Teukolskvl 1 19861 ) . 

With the pressure calculated the next step is to calcu- 
late the column density between the particular density and 
the neutron star surface: P/g^s- Where <?ns is the surface 
gravity given by: 



5ns = 



GM / 
R 2 { 



_ 2GM \~ 1/2 
~ (PR J 



(10) 



For a 10km and 1.4Mq neutron star the surface gravity 
is 2.43x 10 14 cm/s 2 . Finally, the minimum mass fraction re- 
quired for an isotope to rise to form the atmosphere is given 
by the ratio of 1 g/cm 2 to the calculated column density. 

For a mass fraction sufficient to form the envelope we 
are looking for isotopes which have abundances on the or- 
der of parts per million over the entire crust, versus in the 
atmosphere where the isotopes need to have abundances on 
the order of parts per billion or less. In order to calculate the 
minimum mass fractions required for an isotope to float to 
form the envelope we scale according to the column densi- 
ties, we use the typical column density of the envelope about 
4 x 10 9 g/cm 2 and 1 g/cm 2 for the atmosphere. 



4 FREEZE-OUT 

To calculate the neutron star freeze-out we consider three 
different timescales: the cooling time (r c , for example Eq. fT} 
for the modified Urea case), the settling time (t s ), and 
the nuclear reaction timescale (r rxn ). The composition of 
the neutron star depends on the timescales; for the case 
Trxn < r c < r s , we expect the particular species to be in 
nuclear statistical equilibrium (NSE); for r c < r rxn < r s 
the particular nuclear reactions are quenched - the star is 
cooling faster than the reactions can occur; for the case 
t s < Trxn < t c the isotopes can float up faster than the 
reactions bring them to NSE. For the case t s < rf rcoze the 
light isotopes can float all the way up to the top before the 
layer freezes. 

The time when the particular layer of the neutron star 
freezes is determined by comparing the potential energy be- 
tween the ions that compose the crust to their kinetic energy 
i|Shapiro fc Teukolskvl Il986h . 



Potential Energy (Z e) 2 



a T kt' 



(11) 



1.42180 x 10 / >(a;)dyne/cm / 



Kinetic Energy 

where e is the electron charge and a is the ion radius such 
that the product of (4/3)-7ra 3 and the ion number density is 
unity. When F > 180 the layer freezes out, or crystallizes, 
and the light species can no longer float upward. 



4.1 Settling timescale 

A second question is whether the light isotopes can reach 
the surface before the layer freezes or before we look at the 
star. On the other hand the gravitational settling could be so 
efficient that the light isotopes could float up before the layer 
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cools and nuclear reactions quench. For these two reasons an 
esti mate of the settling time is crucial. 

iBrown et al.l (|2002h calculate the sedimentation or set- 
tling timescale for a neutron-star atmosphere, 

7 3.9 7 0.3„1.3 

^ 1n S ^1 ^2 p5 no \ 

S [ 4 1.8 ffl4 2 T 0.3 {MZl _ AlZ2) \ ■ W 

This is an estimate of the time for a nuclide to settle down 
over a pressure scale height — a negative value means that 
the nuclide ascends. This timescale is typically around a 
few years for the envelope and about 10 6 yr for the outer 
crust (p < 10 12 g cm -3 ) for Silicon-28 in a background of 
Iron-56. In particular until the bulk of the Nickel-56 has de- 
cayed the settling time for Silicon-28 is much larger because 
— A1Z2 <C 1; consequently, the nucleons differentiate 
gravitationally after the nuclear reactions effectively cease, 

I.e. T rX n <C 7~s • 

4.2 Reaction timescale 

For the modified Urea case we also calculated the reaction 
rate time scales by making use of the subroutines in torch. 
Each of the rate calculations depends only on the input tem- 
perature and the densities of the various species. As there 
are many different ways to make a specific isotope, e.g 28 Si, 
the rates which lead to the creation of the isotope are added 
together to get the timescale of the reaction rate, r rxn . 

In order to calculate the abundance of the alpha parti- 
cles and the other species we make use of the nuclear statis- 
tical Saha equations as implemented in nse. We assume that 
a particular species freezes out of equilibrium when the reac- 
tion timescale exceeds the cooling timescale. The abundance 
in nuclear statistical equilibrium at the freeze out tempera- 
ture gives an alternative estimate of the final abundance of 
the nuclides. 



5 RESULTS 

In order to determine the expected composition of the 
neutron star atmosphere, in the cases of the modified Urea 
and the thick crust we examined a density of 10 7 g/cm 3 . For 
the thin crust a density of 10 7 g/cm 3 would crystallize before 
any of the isotopes had time to reach the surface, so we ex- 
amined a density of 10 6 g/cm 3 . The results from the nuclear 
reaction network are compared with those of a semi-analytic 
freeze-out calculation in the modified Urea case. Each of the 
three cases are discussed below. 

5.1 Case 1: Modified Urea 

5.1.1 Mass Fractions: Atmosphere 

At a density of 10 7 g/cm 3 the corresponding pressure is: 
1.1 xIO 24 dyne/cm 2 . At this pressure the column density to 
the surface is: 4.4 xl0 9 g/cm 2 . The resulting required min- 
imum mass fraction required for an isotope to be optically 
thick on the surface is: 2.3xl0~ 10 . 

Isotopes with a mass fraction greater than 2.3xlCP 10 
will have a surface density of lg/cm 2 . The lightest elements 
to be optically thick on the surface and have time to reach 
the surface before crystallization of the layer occurs are 




2xl0 9 4xl0 9 
Temperature [K] 

Figure 2. Lightest isotopes for which the mass fraction abun- 
dance would be great enough that the isotope will be optically 
thick on the neutron star surface. All of these isotopes have time 
to reach the surface before the layer crystallizes. These are the 
mass fractions for the density of 10 7 g/cm 3 . The corresponding 
pressure and column density for this neutron star density are 
1.1 X 10 24 dyne/cm 2 and 4.4x 10 9 g/cm 2 , respectively. The hor- 
izontal line indicates the minimum abundance required for an 
isotope to have a surface density of lg/cm 2 . The abundances for 
28 Si, 32 S, and 36 Ar in nuclear statistical equilibrium using nse 
are depicted for comparison by the nearly vertical, heavy lines. 

shown in Figure O where the horizontal line indicates the 
minimum mass fraction required to be optically thick on the 
surface. The lightest elements to rise to the surface which 
are optically thick are: 28 Si, 32 S, 34 S, and 36 Ar. In particular 
the abundance of 28 Si is about 3 x 10 -9 so a layer ~ 10g/cm 3 
of silicon lies on the surface of the star. 

5.1.2 Freeze- Out 

Using the steps outlined in 2] we have calculated the cooling 
(r c ), settling (r s ), and the nuclear reaction (r rxn ) timescales 
for the case of 28 Si. The results of these calculations are dis- 
played in Figure [3] These calculations compare the age of 
the neutron star to the settling, crystallization temperature, 
creation and destruction timescales of 28 Si for two densities: 
10 7 g/cm 3 and 10 12 g/cm 3 . The temperatures at which the 
layers crystallize are 4.8 x 10 7 K and 2.2 x 10 9 K, for the 
densities of 10 7 g/cm 3 and 10 12 g/cm 3 , respectively. The cre- 
ation and destruction rates are d In Xi/dt, where Xi is the 
abundance of the 28 Si isotope. These rates are calculated by 
using the routines in torch to determine the energy release 
per unit mass, these are then multiplied by the abundances 
calculated from the nse code. The abundances were also cal- 
culated using the nuclear Saha equation and are the output 
of the nse routine. The two different methods for calculating 
the relative abundances, the output from torch and using 
nse are displayed in Figure [2] It is clear that the results from 
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Figure 3. Comparing the cooling, settling, and nuclear reaction 
timescales of 28 Si for the neutron star densities of 10 7 g/cm 3 and 
10 12 g/cm 3 . The timescales for the density of 10 12 g/cm 3 are de- 
noted by the heavier line, while the thinner line is for the density 
of 10 7 g/cm 3 . The heavy vertical line denotes the temperature at 
which the layer at 10 12 g/cm 3 crystallizes, the temperature for the 
crystallization of the layer at 10 7 g/cm 3 is cooler than 10 8 ' 5 K and 
so lies to the left of the plot. The layer of 10 12 g/cm 3 crystallizes at 
2.2 X 10 9 K, and the layer at 10 7 g/cm 3 crystallizes at 4.8 X 10 7 K. 
The numbers in the parentheses give the corresponding surface 
density at the two neutron-star number densities. Note that the 
cooling and reaction timescales overlap for the two densities. 



torch do not follow nuclear statistical equilibrium precisely. 

Using Figure [3] we can find the quenching temperature 
for the silicon reactions. The reactions become quenched 
when the reaction rate time becomes longer than the age of 
the neutron star. From the plot the quenching temperature 
for the reactions which destroy silicon are 2.12 x 10 9 K and 
2.15 x 10 9 K for the densities of 10 12 g/cm 3 and 10 7 g/cm 3 , 
respectively. The reactions which create silicon are quenched 
at 3.05 x 10 9 K for both of the densities. These approximately 
give the temperature range over which the abundance of sil- 
icon levels out in Figure [2] The relative abundances of 28 Si 
at the quenching temperatures for the destructive reactions 
are 4.86 x 10" 13 and 2.36 x lO" 10 , for densities of 10 12 g/cm 3 
and 10 7 g/cm 3 , respectively. The relative abundances of sil- 
icon for the creation reactions are: 5.82 x 10~ 9 at a density 
of 10 12 g/cm 3 and 1.84 x 1(T 6 at a density of 10 7 g/cm 3 . 



5.2 Case 2: Thick Crust 

In the case of the thick crust, the density of 10 7 g/cm 3 was 
examined for the isotopes which can rise to the surface. The 
pressure, column density to the surface and the minimum 
mass fraction required to be optically thick are the same 
for the modified Urea case. For the thick crust the lightest 



Figure 4. The lightest isotopes for which the mass fraction abun- 
dances are large enough that the isotopes which can reach the 
surface are optically thick. These are the mass fractions for a 
density of 10 7 g/cm 3 and for a neutron star with a thick crust. 
The horizontal line indicates the minimum abundance required 
for an isotope to have a surface density of lg/cm 2 . 

isotopes which are optically thick and can rise to the surface 
are 50 Cr, 53 Mn, 54 Fe, 55 Fe, and 57 Co, as shown in Figure H 

5.3 Case 3: Thin Crust 

For the thin crust the isotopes which had the possibil- 
ity to be optically thick could not reach the surface be- 
fore the layer crystallized for densities at 10 7 g/cm 3 and 
higher; consequently, we examined the layer at a density 
of 10 6 g/cm 3 . At this density the corresponding pressure is: 
2.3 x 10 22 dyne/cm 2 . The column density to the surface at 
this pressure is 9.6 x 10 7 g/cm 2 , requiring a minimum mass 
fraction of 1.0 x 10 -8 for an isotope to optically thick on 
the surface. The isotopes which are not only optically thick, 
but can also rise to the surface are shown in Figure [S] These 
isotopes include 40 Ca, 50 Cr, 53 Mn, 54 Fe, and 55 Fc. 



6 CONCLUSIONS 

In the case of cooling via the modified Urea process, with- 
out the thermal influence of a crust, our results (e.g. Fig. [3} 
show that silicon has sufficient time to float to the top from 
a density of 10 7 g/cm 3 before the layer freezes or we observe 
it. However, the settling time from 10 12 g/cm 3 is too long 
for light species to float up before the crust freezes; there- 
fore, we can conclude that the atmosphere in this case is 
likely to be composed of silicon but the envelope is likely to 
be composed of iron-group elements that have been chemi- 
cally separated by gravitational settling. Deeper layers are 
unlikely be chemically separated at least by gravity. 

We have used the torch code in order to calculate the 
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Figure 5. The lightest isotopes which can reach the surface 
and are also optically thick. These isotopes are for a density 
of 10 6 g/cm 3 for a neutron star with a thin crust. The cor- 
responding pressure and column density at this density arc 
2.3 X 10 22 dyne/cm 2 and 9.6 X 10 7 g/cm 2 , respectively. The hori- 
zontal line indicates the minimum mass fraction required in order 
for an isotope to be optically think. 



lightest isotopes which would rise to the neutron star sur- 
face and have compared these results to semi-analytic cal- 
culations. Using the torch code we calculated the the mass 
fractions at a density of 10 7 g/cm 3 . We found the lightest iso- 
topes to rise to the surface and be optically thick in the at- 
mosphere are: 2S Si, 30 Si, 31 P, 32 S, 33 S, 34 S, and 36 Ar. On the 
other hand, there is not sufficient time for light elements to 
percolate from high density to form the envelope. However, 
calculations of freeze out at higher densities are required to 
determine the precise composition of the envelope. 

We also did semi-analytic calculations of the freeze-out 
of 28 Si, for the modified Urea case, assuming local thermody- 
namic equilibrium until the cooling rate exceeds the reaction 
rate. We used the rates from torch and the abundances from 
the code nse in order to calculate the rates of nuclear reac- 
tions involving 28 Si; the reactions become quenched when 
the reaction time is longer than the age of the neutron star. 
We found the creation reactions are quenched at a temper- 
ature of 3.05 x 10 9 K for both the densities of 10 7 g/cm 3 and 
10 12 g/cm 3 . The reactions which destroy silicon are quenched 
at 2.12 x 10 9 K and 2.15 x 10 9 K at a density of 10 12 g/cm 3 and 
10 7 g/cm 3 , respectively. The calculated quenching tempera- 
tures of the silicon reactions agree with the results calculated 
using the torch code directly. 

In the case of the thick crust, the atmosphere could 
be formed by 50 Cr rising to the surface from a density of 
10 7 g/cm 3 . For a neutron star with a thin crust and direct 
Urea cooling the layers at 10 7 g/cm 3 and higher crystallize 
before the isotopes could have time to reach the surface. The 



atmosphere in the thin crust case could be formed from 40 Ca 
rising to the surface from a layer at a density of 10 6 g/cm 3 . 

Unless there has been significant accretion either from 
the supernova debris, the interstellar medium or a compan- 
ion, neutron-star atmospheres are unlikely to be composed 
of iron, helium or hydrogen. The type of isotope compos- 
ing the atmosphere depends on the cooling mechanism at 
work in the star: 28 Si for modified Urea, 50 Cr for a thick 
crust, and 40 Ca for a thin crust and direct Urea cooling. 
The formation of these atmospheres can provide additional 
justification, with fallback accretion, for isolated neutron 
stars fit with intermediate atmosphere models. Understand- 
ing how this novel composition of the atmospheres affects 
the neutron-star emission may provide new insights on the 
observed spectra of neutron stars. 
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